Exciton dynamics and Quantumness of energy transfer in the Fenna-Matthews-Olson 

complex 
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We present numerically exact results for the quantum coherent energy transfer in the Fenna- 
Matthews-Olson molecular aggregate under realistic physiological conditions, including vibrational 
fluctuations of the protein and the pigments for an experimentally determined fluctuation spectrum. 
We find coherence times shorter than observed experimentally. Furthermore we determine the energy 
transfer current and quantify its "quantumness" as the distance of the density matrix to the classical 
pointer states for the energy current operator. Most importantly, we find that the energy transfer 
happens through a "Schrodinger-cat" like superposition of energy current pointer states. 
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I. INTRODUCTION 

Recent experiments on the ultrafast exciton dynam- 
ics in photosynthctic biomolccules have brought a long- 
standing question again into the scientific focus whether 
nontrivial quantum coherence effects exist in natural bio- 
logical systems under physiological conditions, and, if so, 
whether they have any functional significance. Photosyn- 
thesis [l[ starts with the harvest of a photon by a pigment 
and the formation of an exciton, followed by its transfer 
to the reaction center, where charge separation via a pri- 
mary electron transfer is initiated. The transfer of exci- 
tations has traditionally been regarded as an incoherent 
hopping between molecular sites @. 

Recently, Engel et al. 0, 0] have reported longlasting 
beating signals in time-resolved optical two-dimensional 
spectra of the Fenna-Matthews-Olson (FMO) 0, Q com- 
plex, which have been interpreted as evidence for quan- 
tum coherent energy transfer via delocalizcd exciton 
states. It transfers the excitonic energy from the chloro- 
some to the reaction center and consists of three identical 
subunits, each with seven bacteriochlorophyll molecular 
sites (the existence of an eighth site is presently under 
investigation). Quantum coherence times of more than 
660 fs at 77 K Q and about 300 fs at physiological tem- 
peratures 3 have been reported. 

Together with recent experiments Q on marine crypto- 
phyte algae, these reports have boosted on-going research 
to answer the question how quantum coherence can pre- 
vail over such long times in a strongly fluctuating physio- 
logical environment of strong vibronic protein modes and 
the surrounding polar solvent. Theoretical modeling of 
the real-time dynamics is notoriously complicated due to 
the large cluster size and strong non-Markovian fluctua- 
tions. It relies on simple models Q of few chromophore 
sites which interact by dipolar couplings and which are 
exposed to fluctuations of the solvent and the protein 
d Q . A calculation of the 2D optical spectrum assuming 
weak coupling to the environment does not reproduce the 
experimental data @ . It also became clear that standard 



Redfield-type approaches fail even for dimers (lfj, lll| . 

In connection with molecular exciton transfer, ex- 
tensions beyond perturbative approaches, such as the 
stochastic Liouville equation Q, extended Lindblad ap- 
proaches [l2j or small polaron approaches [lH , have been 
formulated. Two recent approaches have included parts 
of the standard FMO model. A second-order cumulant 
time-nonlocal quantum master equation found coherence 
times [3] as observed experimentally for the FMO com- 
plex. However, they employed an Ohmic bath in which 
the relaxation time is not uniquely fixed by experiments 
and excluding strong vibrational protein modes. On the 
other hand, a variant of a time-dependent DMRG scheme 
[lU has been applied to a dimer in a realistic FMO en- 
vironment [ill. The solution of the full state-of-the-art 
FMO model with seven localized sites and a physical en- 
vironmental spectrum [l6| including a strongly localized 
Huang-Rhys mode 0] is still missing. In particular, it 
has not been clarified whether a strong vibronic mode 
still allows for the observed quantum coherence times. 

In this paper we close this gap and present numerically 
exact results for the full FMO model with seven sites 
and the spectral density of Ref. [lfjj . determined from 
optical spectra [lTj ■ Upon adopting the iterative real- 
time quasiadiabatic propagator path-integral (QUAPI) 
[HG1 scheme we find coherence times shorter than ob- 
served experimentally. For comparisson we recalculate 
also the dynamics employing the same Ohmic fluctua- 
tion spectrum as Ishizaki and Fleming [l4| . To quantify 
quantum effects for energy transfer, we calculate the en- 
ergy current associated with the transfer dynamics and 
its "quantumness" . It has been argued that quantum 
entanglement could be created in a single excitation sub- 
space during the energy transfer (2fjj . but this form of 
entanglement cannot be used to violate a Bell-inequality 
[2l| and its role for the transfer efficiency is therefore un- 
clear. Wilde et al. [22j used the Leggett-Garg inequality 
to discuss whether quantum effects are relevant but they 
applied a phcnomenological Lindblad approach. Here we 
use a recently developed measure of quantumness based 
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on the Hilbert-Schmidt distance of the density matrix 
to the convex hull of classical states [23fl . taken as the 
"pointer-states" [24| of the energy-current operator. We 
show that energy transfer starts with large current out of 
the initial site and then small currents flow between all 
seven sites. All currents show substantial quantumness. 
Interestingly enough, this implies that the energy trans- 
fer happens through a largely coherent "Schrodinger-cat" 
like superposition of pointer states of the energy current 
operator. 

In the next section we shortly recapitulate the Hamil- 
tonian and the fluctuation spectrum of the site energies 
for the Fenna-Matthews-Olson complex. Then we de- 
termine the population dynamics within the single exci- 
tation subspace. In the fourth section we calculate the 
energy currents between any two sites in the FMO and 
then introduce a measure for the quantumness of these 
currents. Finally we conclude with a short summary of 
our results. 
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FIG. 1. (Color online) Bath spectral densities of the FMO 
complex of C. tep idum: The most realistic form as determined 
up to present [T(j,[l7[ is given by Eq. ([3| (AR, red (upper four) 
lines). The blue dash-dot-dashed (lower) line shows the Drude 
form Gd(oo) used in 



II. FMO MODEL 

The FMO complex is a trimcr consisting of identi- 
cal, weakly interacting monomers each containing 
seven bacteriochlorophylla (BChla) molecular sites which 
transfer excitons. The pigments are embedded in a large 
protein complex. Each of it can be reduced to its two 
lowest electronic levels and their excited states are elec- 
tronically coupled along the complex. Recombination 
is negligibly slow (~ ns) compared to exciton transfer 
times (~ps). Thus, the excitation dynamics is reliably 
described within the 1-exciton subspace. The coupling 
of the seven excited levels gives rise to the Hamiltonian 
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_ (1) 

in units of cm 1 in site representation [131 f° r an FMO 
monomer of C. tepidum. We define the lowest site energy 
of pigment 3 as reference. 

The vibrations of the BChla, the embedding protein 
and the surrounding polar solvent are too complex for a 
pure microscopic description and are thus treated with 
the framework of open quantum systems. They induce 
thermal fluctuations described by harmonic modes Q 
and lead to the total Hamiltonian [lql 



H = H 
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(2) 



3 = 1 « j=l K 

with momenta displacement gj, K , frequency uij, K and 
coupling 1/^ of the environmental vibrations at site j. 



We assume that fluctuations at different sites are identi- 
cal but spatially uncorrelated |25| . 

The key quantity which determines the FMO coher- 
ence properties is the environmental spectral density 
GM = Ej, K (l^ ) | 2 /2^, K )5(w - u jtK ). The most de- 
tailed identification of the bath properties up to date 
has been achieved by Adolphs and Renger They 
used an advanced theory of optical spectra, a genetic 
algorithm, and excitonic couplings from electrostatic cal- 
culations modeling the dielectric protein environment to 
derive the up to present most detailed spectral density 



Here, S = 0.5, S H = 0.22, lo h = 180cm- 1 and 



(3) 



So M = 6.105- 10~ 5 — e V^f + 3.8156 -lO" 5 — e v^F 



with cji = 0.575 cm -1 and u>2 = 2 cm -1 . It in- 
cludes a broad continuous part (see red (upper four) lines 
in Fig. [TJ which for u> — > behaves as super-Ohmic, 
Gar(uj) ~ w 5 , and which describes the protein vibra- 
tions with the Huang- Rhys factor So. It was determined 
from temperature-dependent absorption spectra [l7j . In 
addition, a vibrational mode of the individual pigments 
with the Huang-Rhys factor Sh is included. We have 
added a broadening to the unphysical 5-peak, which is 
justified since the protein is embedded in water as a po- 
lar solvent which gives rise to an additional weak Ohmic 
damping of the protein vibrations. We fix its width to 
7 P = 29 cm -1 which was found for the lowest energy 
peak of protein vibrations in the LH2 complex [2(| . 

In order to recover known results for the transfer 
dynamics of excitations in the Fenna-Matthews-Olson 
(FMO) complex embedded in an Ohmic bath [TJ, we 
additionally consider the spectral density used in Rcf. 
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FIG. 2. (Color online) Time-dependent occupation proba- 
bilities of all seven FMO sites for T = 300 K, 77 K with 
Pn(0) = 1 and pe(0) = 1 for the measured FMO spectrum 
Eq. © 13 . 



FIG. 3. (Color online) Time-dependent occupations of the 
seven FMO sites for T = 300 K and T = 77 K with /9n(0) = 1 
and p66(0) = 1 for the Ohmic spectrum Gu(w). 



[3 given by 



G D (u) 



(4) 



with a Debye cut-off at frequency ui c (see blue dash-dot- 
dashed (lowest) line in Fig. [T|). It includes the reorga- 
nization energy A = 35 cm -1 and the environmental 



timescale 



= 50 fs. 



III. POPULATION DYNAMICS IN THE FMO 

To understand whether the vibrational fluctuations of 
the protein and the pigments allow for the observed co- 
herence times, we apply the numerically exact quasi- 
adiabatic propagator path-integral (QUAPI) to simulate 
the real time exciton dynamics in the FMO complex for 
the realistic bath spectral density ([3]). QUAPI is well es- 
tablished [TH, [l9| and allows us to treat nearly arbitrary 
spectral functions at finite temperatures. We extended 
recently the original scheme to treat multiple environ- 
ments, i.e. separate environments for each chromophorc 
site [H]. 

We calculate the time-dependent populations of the 
FMO pigment sites for the spectral density of Eq. 
We choose T = 300K (physiological temperature) and 
T = 77K (typical experimental temperature). Both the 
pigments BChl I and BChl 6 arc oriented towards the 
baseplate protein and are thus believed to be initially ex- 
cited (entrance sites) [H]]. Thus, we consider two cases, 
i.e., pii(0) = 1 and P66(0) = 1. We focus on the transient 
coherence effects and thus do not include an additional 
sink at the exit site BChl 3. In Fig. [2J we show the 
time-dependent pigment occupation probabilities Pjj(t). 
Identical simulations using smaller width, i.e. 7 P = 5 
cm -1 and 7 P = 1 cm' 1 , for the vibrational mode yield 
identical results for the populations (not shown). For 



P66(0) = I at room temperature, coherent oscillations are 
completely suppressed and at 77 K they last up to about 
250 fs. For pn(0) = 1, coherence is supported longer 
due to the strong electronic coupling between sites I and 
2. At room temperature, it survives for up to about 200 
fs and for 77 K up to 500 fs at most. Thus, coherence 
times are shorter than experimentally observed. We em- 
phasize that coherence features are a transient property 
and, hence, not only the low frequency bath modes, but 
also the discrete vibrational modes are relevant (100 fs 
corresponds to ~ 333 cm -1 ). 

We additionally calculate the time-dependent popula- 
tions of the FMO pigment sites for T = 300K (physio- 
logical temperature) and T = 77K (typical experimen- 
tal temperature) for the FMO with the bath spectrum 
Gn(w) (see Eq. Q). As above we consider two choices 
of initial conditions, i.e., pn(0) = 1 and P6e(0) = 1 and 
focus on the transient coherence effects and thus do not 
include an additional sink at the exit site BChl 3. 

Fig. [3] shows the occupation probabilities of all seven 
sites at T = 300 K versus time. For both initial con- 
ditions, coherent oscillations of the populations of the 
initial site and its neighboring site (j = 2 for pii(0) = 1 
and j = 5 for P66(0) = 1) occur due to the strong elec- 
tronic couplings. They last up to ~ 350 fs. At T = 77 K, 
as shown in Fig. 0(b) and[3Jd), the oscillations persist 
up to ~ 700 fs. Our results coincide with those of Ref. 
[14| . Tiny deviations at short times_arise since the full 



141. The coherence 



FMO trimer is considered in Ref. 
times, however, are not affected. 



IV. QUANTUMNESS OF ENERGY TRANSFER 
THROUGH FMO 



Next, we discuss whether energy transfer dynamics is 
a quantum coherent process. Classical coupled dipolcs 
can show oscillatory energy currents identical to quan- 
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turn mechanical dipolc-coupled systems provided that 
the interaction between the dipoles is weak enough to 
be treated via linear response [29|. In FMO (see Eq. [TJ, 
however, off-diagonal couplings are on the same order of 
magnitude as site energy differences. We use a physi- 
cally well motivated measure of quantumness which can 
be found by first identifying the most classical pure states 
in the problem. In quantum optics, for instances, coher- 
ent states are considered "(most) classical" due to their 
minimal uncertainty and their stability under the free 
evolution of the electro- magnetic field [30j . Next, one 
identifies the set of all classical states with the convex 
hull (i.e. all mixtures) of the pure classical states, as mix- 
ing cannot increase quantumness. This leads in quantum 
optics and for spin-systems to the well established notion 
that classical states are those with a well-defined positive 
P- function [23|, l3ll - l33j . Finally, one defines "quantum- 
ness" as the distances to the closest classical states. 

For FMO we are interested in the "quantumness" of 
the energy current since energy transfer constitutes the 
main function of the FMO complex. A measurement 
of the energy current would yield that for any measure- 
ment result the state of the system would collapse onto 
a pointer state. In general, pointer states are the states 
einselected by the interaction with its environment [24| 
being here the hypothetical measurement apparatus for 
measuring energy current. Pointer states are the most 
classical states in that they arise via decoherence and 
are those that persist when a quantum system interacts 
strongly and for a long time with its environment. In our 
case of a projective von Neumann measurement, they are 
the eigenstates of the quantum mechanical operator cor- 
responding to the energy current measurement. 



A. Energy currents in FMO 

An energy current operator can be derived for a general 
multi-site Hamiltonian [34[ . The expression for an energy 
current operator j(x) can be obtained from a continuity 
equation 



ff(x) + V-j(x)=0, 



(5) 



is a 



34], 



where H(x) is the local energy density. This 
well-established procedure introduced by Hardy 
and which was largely adopted by a large number of 
subsequent papers (see e.g. |35l - |39| .) There is some 
freedom in the definition of local energy density, but for 
a Hamiltonian of the form H = Ti + (1/2) Yli^k 
where Tj and Vi arc kinetic energy and interaction 
energy (of possibly a large number of particles) in region 
i, and Vik the interaction between the particles in regions 
i and those in region k, a natural decomposition is given 
by H = hi = Ti+ (l/2)E fc ^V« = ViH(xi), 

where Vi is a volume of region i. 

For the FMO complex, precise information about the 
the single-excitation sector of the macro-molecule with 



sites i € {1, . . . , N} is known (N = 7). The Hamiltonian 
H = Hfmo is given in tight-binding approximation as 



iV 



H = Y,hik\i){k\ 



(6) 



i,k 



where \i) is a state with the excitation localized on site i 
(see Eq.(l) ). The natural decomposition of that Hamil- 
tonian in terms of local excitations is 



N N 

H = J2h, with hi = -^2(h lk \i)(k\ + h.c.) . 



(7) 



i=l 



k=l 



In order to deduce an energy-current operator from this 
Hamiltonian, consider first a linear chain in ID, where j 
has only one component, denoted by j, with values ji on 
sites i, taken as equidistant with lattice constant a. The 
discretized form of §5$ reads 



dt a 



3i~3i 



(8) 



where j i ' r denote the energy flux in positive x-direction 
on the left of site i and on the right of site i. It is con- 
venient to think of j\ as being evaluated half way 
between sites i — 1 and i (between sites i and i + 1), 
respectively. Current j\ arises from the balance of the 
currents s,_i_>j from site i — 1 to site i and s,_>i_i from 
site i to site i — 1. The latter two currents are always 
defined positive and in the direction indicated by the in- 
dices, whereas ji can be positive or negative, depending 
on which current component dominates. Inserting these 
expressions into ([8[l. we find 



at 



E 

k=i±l 



(■n 



(9) 



This equation is valid for the ID tight-binding model 
with only nearest neighbor couplings. It has the natu- 
ral interpretation that the change of energy density at a 
given site is given by the difference between the sum of 
incoming energy currents and the sum of outgoing energy 
currents. As such, the expression generalizes in straight- 
forward fashion to a general tight-binding model on a 
graph, where each site can be connected to an arbitrary 
number of other sites. Each link can support an energy 
current, and one has thus 



at 



(10) 



The left-hand side of (fT0|) is easily calculated by using 
Heisenberg's equation of motion, 



§ i h i = t -[H,hi}. 



(11) 



A short calculation leads to 
d_ 

dt"' 2h 



T^h, = ■^t } j } j (hi C jhu\k)(l\ - h ik h kl \i)(l\ - h.c.) 



k=£i I 



(12) 
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FIG. 4. (Color online) Selected average energy currents jo, 
(top panels) and corresponding quantumness (lower pan- 
els) for FMO spectral density of Eq. © and T = 300K. Notice 
the rescaling of the currents for better visibility. 



Comparing this expression with ((9]), we identify the di- 
rected energy current 



Si_>& = ( h ik h ki\i) (l\ - h.c.) . (13) 

l 

The (positive or negative) energy current attributed to 
link i — k on the graph is jik = Sj-s-fc — Sk^i- The dimen- 
sion of all these energy currents is energy/time. A cor- 
responding classical energy current phase space function 
can be derived in classical mechanics, if a corresponding 
classical Hamilton function can be established, in a com- 
pletely analogous way by using Hamilton's equation of 
motion to derive a continuity equation for energy trans- 
port. 

We have calculated all energy currents (jik(t)) = 
tx\jikp{t)\ for all sites i and k. The currents out of the 
entrance site towards its strongest coupled neighbor (j'12 
for pn(0) = 1 and j'56 for P66(0) = 1) have by far the 
largest amplitude. All currents show initial oscillations 
which are more pronounced at 77K than at 300K (see 
Fig. 0] and Fig|5]) on times scales comparable to coher- 
ent oscillations discussed above. The currents out of the 
entrance site towards its strongest coupled neighbor (ji2 
for pn(0) = 1 and j 56 for P66(0) = 1) have by far the 
largest amplitude. After this initial phase energy cur- 
rents between all sites are small but finite (as exemplary 
shown by j'13 and jse) and persist up to lOOOfs. Energy 
transfer, thus, starts by quantum coherent population ex- 
change mainly between site 1 and 2 ( or 6 and 5 respec- 
tively) and some transfer to other sites. After this initial 
phase, in which a prelimenary redistribution of energy 
takes place, small currents between all sites will slowly 
bring the system into thermal equilibirum and result in 
an according population / energy distribution. 



FIG. 5. (Color online) Same as FigH but at T = 77K. 



B. Classical states of energy transport 

It has long been appreciated that the most classical 
states corresponding to a given quantum mechanical 
observable are the so-called "pointer states". These 
states are "einselected" by the decoherence process due 
to the interaction of an environment with the system. 
In the case of a measurement, the environment is given 
by the measurement instrument, and the pointer states 
are in this case simply the eigenstates of the operator 
representing the quantum mechanical observable [3, 53 ■ 
It is in the eigenbasis of these pointer states that the 
density matrix becomes diagonal due to the decay of 
the off-diagonal matrix elements. The remaining diago- 
nal matrix elements correspond to the probabilities of 
finding one of the possible outcomes of the measurement. 

As eigenstates of a Hcrmitian operator, pointer states 
are pure states. The most general classical states can 
be obtained by classically mixing the pointer states, 
i.e., one chooses pointer states randomly with a certain 
probability. As this is a purely classical procedure, 
it cannot increase the "quantumness" of the system. 
This definition of "classicality" is well-established in 
quantum optics and leads to a criterion of classical- 
ity based on a well-defined positive-definite P-function 
[il], 51]. It was recently extended to spin-systems f23l l3lj ] . 



C. Quantumness 

Thus, if one wants to decide whether the energy cur- 
rent between two different sites in FMO has to be consid- 
ered "quantum", one has to (?) find the pointer states of 
the corresponding energy current operator, and (ii) check 
whether the state of the system p(t) can be written as a 
convex combination of these pointer states. A more quan- 
titative statement is possible by measuring the distance 
of p(t) to the convex set of classical states @,|3ll,53, i-e., 
here the convex hull of the pointer states. The pointer 
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states are the eigenvectors \v\ k ) of jik- They define the 
relevant pure classical states for the energy transfer along 
link i — k. Their classical mixtures form the convex set of 
all classical states for the energy current on link i—k. Any 
distance measure in Hilbert space is in principle suitable, 
and we use the Hilbert-Schmidt distance for simplicity, 
leading to our definition of quantumncss of the current 

Jik j 

Qik{p)= r , mm \\p-y^Pl\v l ik )(v l ik \\\ (14) 

{Pi\Pi>o,J2iPi =1 } j 

with P| | = (txAA^) 1 / 2 . Note that this measure of 
"quantumness" is completely analogous to and on the 
same level of abstraction as the measure of entanglement 
based on the distance of a state to the convex set of 
separable states (see p. 363 in Ref. [44[), but has the ad- 
vantage of being meaningful even in the single-excitation 
sector, and without artificially separating the system into 
two subsystems. 

A finite value of Qik means that the state cannot be 
written as a convex sum of pointer-states of the energy 
current operator ji k , i.e., there are coherences left in p(t) 
written in the pointer basis. This means essentially, that 
the system is in a Schrodingcr-like cat state of different 
energy-current pointer states at a given time t. 

By definition, Q(p) > 0, and Q(p) — if p is clas- 
sical. An upper bound is given by Q(p) < Q max = 
y/trp 2 — 1/d, where d is the dimension of the Hilbert 
space [lH (Qmax — 0.925 for pure states and d = 7). We 
have calculated the quantumness Qik(t) for all energy 
currents (jik{t)). As discussed above the currents out 
of the entrance site towards its strongest coupled neigh- 
bor (ji2 for pii(0) = 1 and j 56 for p 66 (0) = 1) have by 
far the largest amplitude and all currents show initial 
oscillations which are more pronounced at 77K than at 
300K (see Fig.|4]and Fig. [5]) on time scales comparable to 
coherent oscillations discussed above. These oscillations 
go hand-in-hand with substantial quantumness: Q12 and 
Q13 arc initially of the order of Q ~ 0.6 — 0.8, but drop 
within ~ lOOfs to Q ~ 0.2 — 0.4, with the exception of 
the case (pii(0) = 1, T = 77K), where Q13 drops slowly 
over lOOOfs, much more slowly than the oscillations of 
the energy current, and the rapid decay of coherences in 
the site basis not withstanding. After this initial phase 
energy currents between all sites are small but finite (as 
exemplary shown by and j'36) and persist up to lOOOfs 
and, noteworthy, a substantial amount of quantumness of 



the order Q ~ 0.2 — 0.4 still persists at lOOOfs as well. In 
the last graph (pqq(0) = 1, T = 77K), the quantumness 
even rises again after the initial drop. 

This implies that the superposition of the pointer 
states of the energy current operator remains largely 
coherent even after population dynamics does not show 
coherence anymore. Nature has apparently engineered 
the environment of the FMO complex in such a way that 
it is rather inefficient in decohering superpositions of 
energy currents. Thus even though the actual currents, 
which bring the system into thermal equilibrium, are 
rather small (after the initial oscilatory phase), they 
turn out to have a finite quantumness. This is further 
illustrated by comparing with the quantumness of the 
thermal equilibrium state px = exp(— H-p^o/hsT) at 
T = 77K and T = 300K, which is of the order 0.01 and 
0.002, respectively. Altogether, our data quantitatively 
show the non-classical nature of the energy transfer 
in FMO, and provide a clear physical picture of the 
quantumness: energy transport in FMO happens largely 
through a coherent Schrodinger-cat like superposition of 
pointer states of the energy current operator. 



V. CONCLUSION 

To summarize, we have obtained numerically exact re- 
sults for the real-time exciton dynamics of the FMO com- 
plex in presence of realistic environmental vibronic fluc- 
tuations. We have used the most accurate form of the 
spectral density realized in nature. It includes vibronic 
effects via a strongly localized Huang-Rhys mode. The 
resulting coherence times of the populations are shorter 
than observed experimentally. The energy transfer dy- 
namics is also intrinsically quantum mechanical on the 
time scales of several hundred femtoseconds. This has 
been shown by calculating the quantumness as distance 
to the convex hull of pointer states of the energy cur- 
rent. The energy transport in FMO is to a large extent 
through a coherent Schrodinger-cat like superposition of 
the pointer states. 
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